PROCEEDINGS 

OF SCIENCE 



o 
o 

(N 
-i— > 
o 

O 



Eradication of singularities in the next-to-leading 
order RG evolution for the AS = 1 effective 
Hamiltonian with 3 quark flavours 



Oh! 



> 
OS 



o 



David H. Adams* 

Department of Physics and Astronomy, Seoul National University, Seoul, 151-747, South Korea 
E-mail: [dadams@phya . snu .ac.kr 



Weonjong Lee 

Frontier Physics Research Division and Center for Theoretical Physics, 

Department of Physics and Astronomy, Seoul National University, Seoul, 151-747, South Korea 



E-mail: ^leeSphya . snu .ac.kr 



We consider the renormalization group (RG) evolution for the operators in the AS — 1 effective 
Hamiltonian with 3 active quark flavors, which is needed in the numerical analysis of data sets for 
e'/ £ calculated in lattice QCD. Singularities are present in the original solution of Buras et al. at 
next-to-leading order (NLO). We show how these can be eradicated through a method of analytic 
continuation to obtain the correct finite solution in this case. Furthermore, we trace the origin of 
the singularities to a breakdown of the approach of Buras et al. in the 3 flavour case, and show 
how it can be rectified so that singularitites are absent from the beginning. 
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1. Introduction 

Weak AS = 1 decays of hadrons can be described in the framework of the Standard Model by 
an effective Hamiltonian J^ e ff obtained by integrating out the degrees of freedom associated with 
heavy quarks and gauge bosons ||, ||[ §, §|. It governs interesting aspects of kaon physics; in 
particular the direct CP violation parameter e' /e and the AI = 1/2 rule. The effective Hamiltonian 
has the form 

^ff = ^ECi(M.S,e)Gi(M,S.«) (1.1) 

where d(n,g,e) denotes the Wilson coefficients, the Qi(n,g,e)'s are a basis of renormalized 4- 
fermion operators contributing to the effective Hamiltonian, p, is the scale, g the strong coupling 
and e the EM coupling. We take the basis operators Q\ to be the ones specified by Buras et al. 
in Section 2 of Ref.jjj]]. There are 10 of these and their origins are as follows: <2i-2 are from 
W-exchange, £23-6 are from QCD penguin diagrams while (27- 10 are from electroweak penguin 
diagrams. 

The energy scale in kaon decays is about 500 MeV, so the hadronic matrix elements of ^// 
are dominated by the strong interaction and must be calculated non-perturbatively. The lattice 
approach is the only possibility for doing this at present, and is currently being used; see Ref.'s 
[^, [7[ ||, Q]. The procedure for calculating the matrix elements from the lattice can be summarized as 
follows. First one calculates the matrix elements of the appropriate lattice operators Qj(a). Next 
one matches to the continuum operators at some scale q* 1 /a in a continuum renormalization 
scheme, e.g., Naive Dimensional Regularization (NDR): 

Qi{q*)=z i j{q*a)Q l f t {a) (1.2) 



The matching factors Zij(q*a) can be calculated in perturbation theory [12]; they involve log(^*a) 
so the continuum scale q* must be close to 1 /a to avoid large logs. Finally, one uses RG evolution 
to run the continuum operators £2,(£t) from the scale pt = q* down to pt = m c (the charm quark 
mass) where they can be combined with the known Wilson coefficients 1 at the scale pt = m c to get 
the matrix elements of ^//. 2 

The RG evolution operator required in the final step above has been calculated up to next- 
to-leading order (NLO) by Buras et al. [JTJ| . However, the expression obtained there contains 
singularities in the case where there are 3 active quark flavors. Therefore it cannot be used for 
extracting the matrix elements from lattice QCD simulations with Nf = 2 + 1 dynamical fermion 
flavors. This is a major problem since unquenched lattice simulations with 2+1 sea quark flavors 
are currently underway with a variety of fermion discretizations (see, e.g., [|(]] and references 
therein) and will be addressing kaon physics such as e'/s and the AI = 1/2 rule in the coming 
future. It is therefore imperative to deal with the singularity problem in the NLO expression for 
the RG evolution operator. The full RG evolution operator is known to be singularity-free, so the 
singularity in the NLO expression in the 3 flavor case should be a removable artifact. In this paper 



'The Wilson coefficients are all known up to NLO in peturbation theory 

2 Note that J^° e tt itself is independent of the scale, so we can combine the operators and Wilson coefficients at any 
chosen scale to determine its matrix elements. 
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we review the solution to this problem given recently in Ref.[[Tl]]. The singularity is eradicated by 
a method of analytic continuation to obtain the correct finite NLO expression for the RG operator 
in the 3 flavor case. Furthermore, we trace the origin of the singularities in the work of Buras et 
al. to a breakdown of their approach in the 3 flavor case, and show how it can be rectified so that 
singularities are absent from the beginning. 



2. Review of RG evolution and the singularity problem at NLO 

The running of the <2,-(ju)'s and C,(/i)'s is governed by the RG evolution equation resulting 
from j^J^eff = 0. Combining the operators and Wilson coefficients into vectors Q and C respec- 
tively, the evolution is given by 10 x 10 matrices acting on these vectors. The evolution of <2(a0 is 
obviously inverse to that of C(jx) so it suffices to determine the latter, which is what we do in the 
following. The running of the EM coupling e is negligible over the range of scales that are relevant 
for the present considerations, so we treat it as constant in the following, as was done in Ref-tQ]. 

The RG equation for C(ju) is 

^ + P( g ,e)-^]c = f(g,e)C (2.1) 

where y(g,e) is the 10 x 10 anomalous dimension matrix (given below) and fi(g,e) is the beta- 
function, given by 

« s ' c > = -*T&-ft(W-' i ''(iW + - <2 " 2) 

with 

? 38 8 d 

A, = ll--/, ft = 102-y/, ft e = -|( B +|) (2.3) 

where u and d denote the number of active u— and d— type flavors, respectively, and f = u + d 
(= Nf) is the total number of active flavors. 

From the RG equation ( |2~l| ) the running of C(ju) is found to be given by 

C(mi) = U(mi,m2)C(mQ,) (2.4) 

where the evolution matrix is 

U(m h m 2 ) = T.exp ( [ 8m dg'fj^fl) (2.5) 

with the dependence g = g{m) specified by ra^| = j3(g,e). Here T g denotes g— ordering; it is 
required since generally [y{gi),y(g<i)] ^ for gi ^ g 2 . 

To evaluate U{m\,m2) we need to know the anomalous dimension matrix y(g,e). It is de- 
termined from the renormalization constant matrix relating the bare and renormalized operators: 
Of' =z y Gu)e;(M)and 
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This depends on g and e through a s = and a = and has been calculated perturbatively up to 2 
loops in the NDR scheme, whereby the matrices in the following expansions have been determined 
(see ^ and the references therein): 



where the pure QCD part is 



Ys(g 2 



a 



+ 4^ rfe) + 0(a 



« s (0) , 
4^ + 



(4k) 



and the leading order QED correction in ( |2.7| ) is specified by 

,2\ 



For later use we note that 

y(g,e) 4k 



P(g,e) 



Pog 3 

J°) 



i-. (0) , a s (1) 

^{g)=Ye +-^Tse 1 



(2.7) 



(2.8) 



(2.9) 



+ « ( ri 0) + - £yr - fttf u ') + o(«, 2 ) ) + o(« 



„(1) ft JO) ft, JO) 



= ite +0(s,+0(K) 

The expansion of the evolution operator in the EM coupling takes the form 

U imx^m-i) = UJmi,m2) + —R{m\,m2) + 0((X 2 ) 

4% 



where the pure QCD evolution is 

U s (mi,m 2 ) = T g exp I / dg 



■77 „n 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



and the leading additional contribution to the QCD evolution in the presence of EM interactions in 
( J2. 12[ ) is given by (see 



R(mi,m2) 



g{m ) d , U s (m u m')T T (g')U s (m',m 2 ) 
g(m 2 ) ft(g') 



(2.14) 



where g' = g'(m'). The EM contribution to the beta-function has been ignored in ( 2. 13 )— ( 2.14 ): 
p s (g) = j3(g,0) so the expressions are valid when the pi e term in (|Z2|) is dropped, which is a 
justifiable approximation made in Ref.[[T]]. However, the generalization of the NLO expressions for 
the evolution matrix to the case where the j3i e term is not dropped is straightforward: In light of 



YsP — % yi°^ hi the relevant expressions [pT|]. 



(2. 10) it can be obtained simply by replacing y 4 

In the remainder of this paper we restrict our attention to the pure QCD evolution U s (m\,m2) 
which is where the aforementioned singularity problem arises at NLO in the 3 flavor case. Once 
the singularity is eradicated, the new finite expression needs to be used in the NLO evaluation of 
(2.14) for R(m\,m2). We omit that part here, referring to Ref.[ll] for the details and result. 
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We begin by recalling the leading logarithmic approximation to U s (m\,m2): it is obtained by 



keeping only the leading contribution to in the integrand in (2.13), which is shown in (2. 11). 



It is useful to diagonalize js 



.(0)r . 



yffl = V 1 Js°^ T V diagonal matrix 
then the leading logarithmic approximation is found to be 



TT (o)( s ^ a,(m 2 ) ^ 0) /2fi> t 
Us '(m u m 2 ) = V[ — ; — rl V 



.a s {mi)J 

The NLO contribution to the evolution matrix can now be found starting from the ansatz 



U. 



(mi,m 2 ) = ( 



1 + 



aj/ni) 
An 



j)U s i0 \m u m 2 )(l 



OC s (m 2 ) 
An 



(2.15) 



(2.16) 



(2.17) 



and using the RG equation for U s (m\,m2) (given by (gJJ) with fi(g,e) and y(g,e) replaced by f5 s (g) 
and Ys(g)> respectively) to derive an equation for /. From solving this equation at lowest order in 
a s Buras et al. find ^ U 



where 



j = VSV~ l 



.(o) Pi 



2/3 2 2Po(l + a t - aj) 



with 



a.i = Vth diagonal element of Yp' /2/3q 



(2.18) 



(2.19) 



(2.20) 



From (2.19) we see that the resulting NLO expression for the evolution matrix has a singularity if 
1 +a, — aj = 0. As mentioned in [JIJ], this happens for (z',j) = (8,7) when there are 3 active quark 
flavors (f = 3). 



3. Solution of the singularity problem via analytic continuation 



The expression ( |2.17| ) can be written up to NLO as 



U^ LU (m l ,m 2 ) = U^{m u m 2 ) + ^-VA{m l ,m 2 )V 

An 



1 



-l 



where 



VA(m\,m 2 )V 1 = a s (mi)JUs\mi,m2) — a s (m 2 )Us°\mY,m 2 )J . 
Inserting the expression (2.16) for uf\m\,m 2 ) leads to 



Aij(mi,m 2 ) =Sj 



/a 2 \ a i /cc 2 \ ai 

05i — -0C 2 [ — 



(3.1) 



(3.2) 



(3.3) 
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where at = tt s (mk) for k = 1,2. In the singular case (f = 3, i = 8, j = 7) we now regularize by 
replacing 



a,j + e 



(3.4) 



in ( |2. 19| ). The regularized quantity then becomes 5y 
making the same regularization ( f3~4l ) there, we find 



2&)s- Inserting this into (3.3), and 



Aij(m u m 2 ) 



\2p J e 



/a 2 \ a i+ E 



V2pV £ tt2 voci 



a 2 



e log(^)+0(e 2 ) 



(3.5) 



Thus a finite expression is obtained in the limit where the regularization is lifted. Substituting this 



into ( |3.1[ ) we get a finite NLO expression for the evolution operator. This solution of the singularity 
problem clearly amounts to analytic continuation of the evolution operator (regarded as a function 
of the matrix y^). 



4. Understanding the "singularity" from first principles 

In the following it is convenient to use the notation U s (gi,g2) rather than U s (m\,m2); this is 
justified since the dependence on m 1 ,m 2 enters exclusively through g{mi),g(m 2 ). 

By general arguments the full evolution operator (for pure QCD) can be written as 

Us(gug2)=(i+^j(gi))u^\g 1 ,g 2 ){i+^j(g 2 )y 1 (4.D 

Using the RG equation 

d Y T (s) 

—U s (g,g Q ) = ff^U s (g,g ) (4.2) 
dg PA8) 

to derive an equation for J(g) at leading order, one finds different equations in the "singular" case 
1 + cij — aj = and "non-singular" case JTT|]. 

In the non-singular case, taking J(g) =J + 0(g) leads to a consistent solution for the constant 
matrix /. This is the solution found by Buras et al. that we reviewed in §2. On the other hand, in the 
singular case, the ansatz J(g) =J + 0(g) does not admit a consistent solution. In fact, it turns out 
that Sij(g) = (V~ l J(g)V)ij must diverge for g —> 0. Specifically, in the singular case 1 + at — aj = 
we found in Ref.fllTl] that the leading order equation for Sjj(g) is as follows: 

Pog-fSij(g)=-Gij (4.3) 

The solution is 

Sij(g) = -^og(g)+Cij = -^log(ofc) + 4 (4.4) 
Pa -ipo 
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where c,y is an undetermined integration constant and c- y - = c, ; - — 2^1og(47i). Thus Sy(g) diverges 
for g — > as claimed. Note however that g 2 Sij(g) vanishes for g — > 0, implying that g 2 J(g) vanishes 
in this limit as it should in order for (4.1) to reduce to Uj, ■ (mi,m2) in the small g limit. We also 



note that Sn (g) is actually a function of g 2 (or oQ as it should be. Finally, our previous expression 



(3.5) for Aij(mi,mz) is readily reproduced from (4.4) (see [11]), so the first principles solution 
of the singularity problem presented here agrees as it should with the finite NLO result for the 
evolution operator obtained via analytic continuation in the previous section. 

5. Conclusions 

We have eradicated the singularities in the original solution of Buras et al. to get a finite 
expression for the RG evolution matrix at NLO in the 3 flavor case. This is essential for being able 
to evaluate the matrix elements of the AS = 1 effective Hamiltonian from the lattice with 2+1 sea 
quark flavors. 

The breakdown of the ansatz J(g) = J + 0(g) in Buras et al.'s approach, and the rectifica- 
tion discussed here, illustrate some general subtleties to bear in mind when evaluating evolution 
operators in general. 
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